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Graphene SU(4) quantum Hall symmetry is extended to SO(8), permitting analytical solutions for graphene 
in a magnetic field that break SU(4) spontaneously. We recover standard graphene SU(4) physics as one limit, 
but find new phases and new properties that may be relevant for understanding the ground state. The graphene 
SO(8) symmetry is found to be isomorphic to one that occurs extensively in nuclear structure physics, and 
very similar to one that describes high-temperature superconductors, suggesting deep mathematical connections 
among these physically-different fermionic systems. 
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Graphene in a strong magnetic field has approximate SU(4) 
symmetry M, which permits examining explicit symme¬ 
try breaking by small terms in the Hamiltonian. However, 
the ground state is strongly insulating, with rapid divergence 
of longitudinal resistance at a critical magnetic field |@]. 
The dependence of on sample impurities suggests that the 
resistance is an intrinsic property of an emergent state dif¬ 
fering qualitatively from perturbed SU(4) solutions (sponta¬ 
neous symmetry breaking) SU(4) symmetry can suggest 
the form of possible emergent states but cannot describe them 
quantitatively. Numerical simulations find various possible 
ground states having similar energies but differing structure. 
Thus the nature of the insulating ground state remains elusive. 

Here we show that SU(4) symmetry can be extended to an 
SO(8) symmetry that recovers graphene SU(4) physics, but 
that implies new low-energy modes that transcend SU(4) sym¬ 
metry and for which solutions may be obtained analytically. 
As a first application we revisit the nature of the ground state 
for undoped monolayer graphene in a magnetic field. 

Good reviews of graphene physics are available film [II; 
we recall here only features relevant for the present discus¬ 
sion. Graphene is bipartite with sublattices A and B; the quan¬ 
tity specifying whether an electron is on the A or B sublattice 
is termed the sublattice pseudospin. The dispersion computed 
in tight-binding approximation | IT iQ indicates two inequiv¬ 
alent sets of points in the Brillouin zone, labeled K and K'. 
The two-fold K degree of freedom is termed valley isospin. 
Near these K-points the dispersion is linear, leading to Dirac 
cones. For undoped graphene the Fermi surface lies at the 
apex of the cones, where the level density vanishes and the ef¬ 
fective electronic mass tends to zero. Hence, near the K points 
low-energy electrons for undoped graphene in zero magnetic 
field obey a massless Dirac equation and behave mathemati¬ 
cally as massless chiral fermions, with chirality related to pro¬ 
jection of the sublattice pseudospin. 

In a magnetic field the massless Dirac equation may be 
solved with an appropriate vector potential and the resulting 
Landau levels (LL) are labeled by integers. The « = 0 level is 
unusual in that it is half filled in the ground state of undoped 


graphene, leading to the anomalous counting observed in the 
graphene quantum Hall effect lllllll. For low-energy ex¬ 
citations in each valley (K or K'), inter-valley tunneling may 
be ignored and the electrons in the valley reside entirely on 
either the A or B sublattice, implying that for the n — 0 LL 
valley isospin and sublattice pseudospin are equivalent labels. 
We shall be concerned primarily with this « = 0 Landau level, 
which has, in addition to the Landau orbital degeneracy, a 4- 
fold degeneracy corresponding to spin and valley isospin. 

The largest energy scales are the LL separation and 
Coulomb energy. For neutral graphene the LL separation is 
approximately three times larger than the Coulomb energy, 
which is in turn much larger than other interactions. Hence, 
we shall ignore inter-LL excitations and consider only a single 
n = Q LL. Justification and caveats for this approximation are 
discussed in Ref. iH]. We adopt a Hamiltonian Hi 

H = Ho-H^+\ Y, [gzi'^i+g±ii'^i + Tj.T;!')] Sin-Fj) (1) 
ifj 

where the Pauli matrices Za operate on valley isospin, the 
Pauli matrices Oa operate on electronic spin, and g± 
are coupling constants, fis is the Bohr magneton, and the 
spin z direction is assumed aligned with the magnetic field. 
The three terms in Eq. ([1]) represent the valley-independent 
Coulomb interaction, the Zeeman energy, and the short-range 
valley-dependent interactions, respectively. 

The four internal states representing possible combinations 
of the projection of the spin a and the projection of the valley 
isospin T are displayed in Fig.[T] Symmetries of the Hamilto¬ 
nian O may be investigated by introducing the 15 operators 

Sa=5I {(y'\C7a\<y)cl^,^^CT:ami, (2a) 

nil, taa' 

Ta = 51 51 (2b) 

ini, atr' 

= SL L (2c) 

ink aa'x 

^ali = jY L (2d) 

ink aa'zz' 
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FIG. 1: Isospin-spin quantum numbers and basis vectors. 


where c^(c) create (annihilate) fermions, a = (x,y,z), j3 = 
(x,y), T and a are defined in Fig.[T] and nij^ labels orbitally- 
degenerate LL states. Physically, is total spin, Ta is total 
valley isospin, Na is a Neel vector measuring the difference in 
spins on the A and B sublattices, and the Hap couple spin and 
valley isospin. Under commutation the operators (|2]i close an 
SUW) algebra that commutes with the Coulomb interaction 
If terms 2 and 3 in Eq. ([T]i are small compared with 
the first, the Hamiltonian has approximate SU(4) invariance. 
Explicit breaking of SU(4) depends on the values of g^ and j 
Eour symmetry-breaking patterns have been discussed |l5l 

Eor a 2-N dimensional fermionic space the most general 
bilinear products cjcj of creation-annihilation operators and 
their hermitian conjugates generate an SU(2A) Lie algebra 
under commutation. Adding the most general pair operators 
cjcy and acj extends SU(2A) to SO(4A) ifTslll^ . The ex¬ 
tended symmetry permits defining a (collective) subspace of 
the full Hilbert space spanned by products of pair creation op¬ 
erators acting on the pair vacuum. An effective Hamiltonian 
constructed from a polynomial in the Casimir invariants of 
all groups in the subgroup chains of SO(4Al) will then repre¬ 
sent the most general Hamiltonian for the collective subspace, 
and will be diagonal in the subspace basis for each dynam¬ 
ical symmetry chain. Thus, the manybody problem can be 
solved exactly in the symmetry limits defined by each sub¬ 
group chain Jy^, and analytically in coherent-state approx¬ 
imation 1 II 8 U 22 I 1 otherwise. This approach has been applied 
extensively to strongly-correlated fermions in various fields; 
for representative examples see |I3, 231-261. 

Eor graphene we assume a single n = 0 LL with creation 
operators and hermitian conjugates Czamt- Degeneracy 
of the LL is denoted by 2Q.. Accounting for 4-fold spin-valley 
degeneracy, 2Q. = 4(20^.) = ABSj{h je), where 20^. is the LL 
orbital degeneracy, B is magnetic field, and S sample size. The 
fractional occupation of the LL is / = «/2n, where n is elec¬ 
tron number, and ths filling factor is v = 4(/— j). 

Now we add to the 15 SU(4) generators of Eq. (I2]i the 
charge operator So = j{n — Q.), the 6 pairing operators 5^ and 
d 1 {ji = 0 , ±1 , ± 2 ), and their 6 hermitian conjugates, with 
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FIG. 2: SO(8) subgroup chains with group generators. Dashed 
boundaries and darker shading indicate subgroups defining the SU(4) 
quantum Hall model. We see that SO(8) subsumes graphene SU(4) 
but has a richer structure with additional subgroup chains. 
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FIG. 3: Configurations created by the operators of Eqs. @ and l|4j 
operating on the pair vacuum |0). Location of the dots (K or K' site) 
indicates the valley isospin; arrows indicate the spin polarization. 


where creates a pair of electrons, one in the a = (ti , ffi) 
level and one in the b = (T 2 ,(T 2 ) level, with the total of 
each pair coupled to zero term by term: A^^ = ' 

We also introduce for later use the linear combinations 

qI = \{S^±dI). (4) 

The 28 operators {§„, Ta, Aa, Ha^;, Hay, 5o, 5, 
generate an SO(8) algebra with a graphene SU(4) subalgebra. 
The full structure for SO(8) subgroup chains is given in Fig.|2] 
Pair configurations created by generators of Eqs. Q and (|4]i 
aerating on the pair vacuum are given in Lig. [3] Kharitonov 
pII has classified collective modes for the n — 0 LL using 
pairs similar physically to ours: creates spin-singlet 

charge density waves (CDW), creates ferromagnetic 

(LM) states, and Qj-lO) creates antiferromagnetic (AL) states. 
These states may be classified according to (§-), which mea¬ 
sures net spin and characterizes FM order, (Tf), which mea¬ 
sures the sublattice charge difference and characterizes CDW 
order, and (Nf), which measures the sublattice spin difference 
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and characterizes AF order. Thus pairs in Eqs. ©-(IHi define 
modes already discussed ||3,S], but now SO(8) symmetry per¬ 
mits analytical solutions for corresponding collective modes. 

Broken-symmetry states based on graphene SU(4) have 
been expressed in terms of the pair wavefunction || 2 l 

IV^)=n E ( 5 ) 

xa-c'a' 

where the product is over the LL orbital degeneracy label mj, 
and the sum is over spin and valley labels. Now let us consider 
SO( 8 ) pairs. All states of an irreducible representation may 
be constructed by successive application of raising and lower¬ 
ing operators to a highest-weight (HW) state (Cartan-Dynkin 
method) Hi]. Let u denote the number of broken pairs. For 
M = 0 states at half filling the pair number isA = jTl = 2k-\-\ 
and the U(4) representation is (^,^,0,0). We choose the HW 
state as the pair state that results from placing one electron in 
the a = \ and one in the a = 2 basis states (see Fig.[T]i, 

in") = m (''l/1"> = (« 

tUk 

where the sum runs over the N states in the LL labeled by 
the OTj. = {—k, —k +l,...,k— l,k) orbital quantum number. 
Writing the sum over in Eq. (| 6 ]l out explicitly and invoking 
antisymmetry eliminates most terms and leaves 

1 / \N mii=+k 

io)= n ^L,4jo). (7) 

mk=-k 

Thus the SO(8) u = 0 HW state is equivalent to a product of 
pairs, one for each nijt^ in the LL. 

Other states can be constructed by applying successively to 
|HW) ladder operators that are functions of the generators G = 
{^a,Ta,Na,Tl„i^} of Eq. For for an arbitrary state |i/) in 
the weight space ji//) = L(G)|HW), where F{G) is specified 
by the Cartan-Dynkin procedure. For example, applying the 
isospin lowering operator T = F{G) = 5 ( 7 ^ — iTy) gives 

I vr) = r_|HW) = n |0)- (8) 

mk 


Likewise, all other states of the u — 0 representation can be 
constructed by using successive applications of raising and 
lowering operators fashioned from the generators of Eqs. 0, 
and they will take the product of sums form (l8]l, just as for 
Eq. 0. Hence the SO(8) D SU W wmmetry introduced here 
recovers existing understanding ISl 1^ of states expected from 
spontaneous breaking of SU(4) by short-range correlations. 

Figure in contains 7 subgroup chains. Each defines a dy¬ 
namical symmetry realized for specific choices of the SO(8) 
Hamiltonian parameters, and yields an exact manybody solu¬ 
tion using standard techniques. We shall deal with these ex¬ 
act solutions in future papers. Here, we interpret the states 
implied by Fig. |2] using coherent state (CS) approximations 
1 221] . The full coherent state solution will be presented in a 


later paper, but we illustrate here for subgroup chains contain¬ 
ing SO(5) in Fig.|2] corresponding to solutions that are linear 
combinations of the symmetry-limit solutions for the SO(8) D 
SU(4) D SO(5) D SU(2), the SO(8) D SO(5) x SU(2)p D 
SO(5) D SU(2), and the SO(8) D SO(7) 3 SO(5) D SU(2) 
dynamical symmetries. The required group theory is already 
known lll7ll22ll27l43l|] . so we can (with suitable change of no¬ 
tation and basis) simply transcribe many equations and rein¬ 
terpret them in terms of graphene physics. Details follow in a 
later paper but we show here results central to this paper. 

Energy surfaces: The energy Eg{n,l5) depends only on n 
and a single order parameter j3. For the symmetry groups g 

E,{n,l5)=N, [A,I5^ +B,{n)p^ + C,{n)+D,{n,l5)] , (9) 

where the group-dependent coefficients are tabulated in Ref. 

. This energy surface results from minimizing 

(//) ~ Go (Csu( 2 )p ) + ^2 (Csu(4)), (10) 

where () is taken in the CS, Cg denotes the quadratic Casimir 
invariant of g, and b 2 and Go are coupling strengths. 

Order Parameter: The order parameter j3 distinguishes the 
phases associated with the subgroup chains in Fig.|2]that con¬ 
tain SO(5). j3 measures AF, since it is related to the AF order 
parameter {Nf) by {Nf) = 2Q\b2\{f - 4Y j3. 

Fluctuations: Coherent states violate translational, rota¬ 
tional, and gauge invariance. However, for realistic fields 
and sample sizes these violations are negligible, yielding a 
Ginzburg-Landau type theory with microscopic pedigree. 

Wavefunctions: Closed forms are given for the SO(8) CS 
wavefunctions in Ref. s. Evaluating these expressions in 
the SO(5) X SU(2) and SU(4) limits, respectively, gives 

ISO 5 X SU 2 ) ~ 1 SU 4 ) ~ (elf |o). ( 11 ) 

Thus the SU(4) state is a superposition of Q± pairs, each with 
vanishing (S^) and {Tfj but finite AF order {Nf). Conversely, 
the SO(5) X SU(2) state is a superposition of S pairs, each 
with vanishing (S^), {Tf), and {Nfj. The critical SO(7) state 
is realized in the transition from SO(5) x SU(2) to SU(4) and 
represents a complex mixture of these wavefunctions. 

Energies for the SO(5) x SU(2), SO(7), and SU(4) limits 
are shown for several values of / = n/2£l in Fig. |4ta)-(c). 
The solutions are distinguished by the AF order parameter 
j3 at the minimum, which is zero for SO(5) x SU(2), non¬ 
zero for SU(4), and indeterminate in the SO(7) critical dy¬ 
namical symmetry that interpolates between SO(5) x SU(2) 
and SU(4) through fluctuations in AF order. For undoped 
graphene the ground state corresponds to the / = 0.5 curves. 
These are shown in Fig.|4jd)-(f) for the three symmetry lim¬ 
its, along with a physical interpretation of the states in terms 
of the wavefunctions (fTTl i. Thus the SO(8) dynamical sym¬ 
metry limits illustrated in Fig. |4ld)-(f) represent a rich set of 
collective states that can be distinguished by the expectation 
value and fluctuations associated with the order parameter j3. 
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FIG. 4: (a)-(c) Coherent state energy surfaces as a function of AF order jS for three SO(8) dynamical symmetry limits. Curves labeled by 
fractional occupation / = n/2Sl (particles or holes, since the SO(8) theory is particle-hole symmetric), (d)-(f) Ground-state energy for three 
symmetry limits. Diagrams indicate the wavefunctions suggested by Eqs. ( lilt , (g) Energy as a function of coupling strength ratio q = Go/bi 
for / = 0.5. Solid green curves (q ^ 0 — 0.5) indicate SO(5) x SU(2) symmetry, solid blue curves for (jr > 1.5 indicate SU(4) symmetry. 
Dashed red curves for g ~ 1 correspond to the critical SO(7) symmetry mediating the quantum phase transition from SO(5) x SU(2) to SU(4). 
(h) Energy surfaces for different occupation fractions / at fixed Gq and b 2 (q = 2.5). Solid green curves for / ~ 0 — 0.2 indicate SO(5) x SU(2) 
symmetry. Solid blue curves for / ~ 0.3 — 0.5 indicate SU(4) symmetry. Curves near / ~ 0.25 (dashed red) correspond to SO(7) symmetry 
mediating the SO(5) x SU(2) SU(4) quantum phase transition. 


Quantum phase transitions between symmetry limits may 
be studied by varying coupling. We rewrite Eq. ( fTOl i in 
terms of a parameter q = bilGo favoring SO(5) x SU(2) 
when ^ << 1, SU(4) when ^ >> 1, and SO(7) when ^ ~ 1 
[(Csu( 2 )p) + (Csu( 4 )) ~ (Cso( 7 )), implying SO(7) symmetry if 
Go b 2 \. Variation of ground state energy with q is shown 
in Fig.Htg). Alternatively, at fixed q phase transitions may be 
initiated by changing particle occupancy. Figure|4jh) displays 
a transition from SO(5) x SU(2) with j5 —0, through a crit¬ 
ical SO(7) symmetry with energy highly degenerate in j3, to 
SU(4) with j3 ^ 0, as / is changed at constant q. 

Thus SO(8) describes analytically a host of broken-SU(4) 
candidates for the states in graphene being unraveled in mod¬ 
ern experiments These solutions provide a spec¬ 

trum of excited states as well as ground states. We shall not 
discuss that here, except to note that all ground state solutions 
have a gap to electronic and collective excitations. The gen¬ 
eral theory to be discussed in forthcoming papers can accom¬ 
modate FM, CDW, and AF states, but for dynamical symme¬ 
tries containing SO(5) all solutions may be classified by a sin¬ 
gle parameter j3 that measures AF order: SU(4) states have fi¬ 
nite p and AF order, but no CDW or FM order, SO(5) x SU(2) 
states have j3 = 0 and no AF, CDW, or FM order, and SO(7) 
states define a critical dynamical symmetry that interpolates 
between SU(4) and SO(5) x SU(2) with no AF order but large 
AF fluctuations, and with no CDW or FM order. We have ne¬ 
glected Zeeman coupling here but it is expected to be small 
for the n = 0 LL ||3l, primarily leading to AF canting J^]. 

Transport properties are not manifest in the algebraic the¬ 
ory but the CS approximation is equivalent to symmetry- 


constrained Hartree-Fock-Boboliubov (HFB) theory |l22 , 


suggesting that SO(8) theory can be mapped onto Hartree- 
Fock (HF) transport calculations. HF calculations for arm¬ 
chair nanoribbons found that AF and CDW states similar to 
ours have no edge currents d. We speculate that our AF 
states also are insulating and thus strong candidates for the 
high-field ground state, but confirmation requires more work. 

Solutions depend on Go and b 2 in Eq. (fTOl i. which define 
effective interactions in the truncated space [highly renormal¬ 
ized relative to parameters in Eq. O]. They may be fixed by 
systematic comparison with data, enabling a robust prediction 
for the nature of the ground and other low-energy states. We 
expect modest impurity levels to modify the effective interac¬ 
tion parameters but leave dynamical symmetries intact. 

The present ideas are similar to ones found in nuclear 
physics It^ and high-Tc superconductors (SC) with 

all three cases exhibiting SO(8) D SU(4) SO(6)] symme¬ 

try and a critical dynamical symmetry generalizing a quantum 
critical point to a quantum critical phase linking other phases 
through quantum fluctuations. In graphene and nuclear SO(8) 
the critical symmetry is SO(7), which interpolates between 
SO(5) X SU(2) and SU(4) il; in SC it is SO(5), which in¬ 
terpolates between SU(2) SC and SO(4) AF Mott insulators 
ll25[ 12611 . These similarities may have implications for cross- 
disciplinary understanding of quantum phase transitions. 

In summary, we have introduced an SO(8) model of mono- 
layer graphene in a magnetic held that recovers SU(4) quan¬ 
tum Hall physics but implies new collective modes transcend¬ 
ing explicitly-broken SU(4) that are leading candidates for the 
high-field ground state. Graphene SO(8) is isomorphic to a 
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symmetry describing many complex nuclei and very similar 
to one describing high-rc superconductors, suggesting a deep 
mathematical connection among these phenomena. 
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